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THERMAL CHARACTERISTICS OF A 
STORAGE VESSEL ON THE MOON 


SUMMARY 

A theoretical analysis has been performed to study the thermal character- 
istics of a liquid hydrogen storage vessel on the moon. Investigations were con- 
ducted concerning the influences on temperature and propellant venting by such 
factors as insulation thermal conductivity, nondiffuse lunar radiation, the shad- 
ow cast by the vessel on the lunar surface, shape and size of the vessel, number 
of isothermal elements into which the vessel surface is divided, and the solar 
absorptance and infrared emittance ratio of the vessel surface. 

No calculations were made to estimate the heat transfer through struc- 
tural penetrations. Neither was an effort made to optimize the weight of the in- 
sulation and the evaporated propellant or to minimize the boil -off by use of tech- 
niques such as shadow shields. 

The factors found to have the most influence were: thermal conductivity 
of insulation, size and shape of vessel, solar absorptance and infrared emit- 
tance ratio, and shadow cast by the vessel on the lunar surface. 

INTRODUCTION 


A computer analysis has been performed to study the thermal character- 
istics of objects on the lunar surface. Because of the probability of storing 
liquid hydrogen on the moon in the not too distant future, a vessel containing 
this liquid was chosen as a typical example. Investigations were conducted 
concerning the influences on temperature and propellant venting by such factors 
as insulation thermal conductivity, nondiffuse lunar radiation, the shadow cast 
on the moon by the vessel, shape and size of the vessel, number of isothermal 
elements on the vessel surface and the solar absorptance and infrared emittance 
ratio ( a / e ) of the vessel surface. 

No calculations were made to estimate the heat transfer through struc- 
tural penetrations. The structural details of the vessel were not defined well 
enough to make a realistic analysis. Neither was a special effort made to opti- 
mize the weight of the insulation and the evaporated propellant or to minimize 
the boil-off by use of techniques such as shadow shields. The computer program 
can handle three vessel shapes: a sphere, a cylinder with hemispherical ends 
and a cylinder with flat ends. Calculations can be performed for any location 
of the vessel on the lunar surface. 



ASSUMPTIONS 

The following conditions were imposed on the analysis: 

1 . The moon is a flat circular area of radius r. 

2. The intensity equator’ coincides with the selenographic equator. 

3. The cryogenic liquid obeys the laws of classical thermodynamics. 

No reduced gravity or nonequilibrium effects are included. 

4. The storage container radiates to the lunar surface and to space and 
receives radiation from the lunar surface and from the sun. No other sources 
or sinks are present. 

5. The storage container is covered with a high performance, multi- 
foil insulation located between metal surfaces. The design of the container is 
only crudely specified. The basic shapes are shown in Figure i. 

6. The ullage pressure at which venting begins is one atmosphere. 
Venting occurs at constant temperature and pressure. The initial temperature 
and pressure are 20. 4°K and one atmosphere, respectively. 

7. An initial ullage of 10 percent, by volume, was assumed. 

THERMAL ANALYSIS 


The temperature of objects on the moon will be nonisothermal, especially 
large surfaces which are only partly illuminated. To minimize this, the in- 
sulation on the storage vessel has been divided into elements. By increasing 
the number of elements each one approaches an isothermal condition at any 
instant of time because its size diminishes. This increases the accuracy of 
the results as well as the complexity of the problem. 

The temperature of an element at any time is given by 


* Line containing subsolar point and subearth point. 
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dT. 

] 

i dt 
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+ D^fS + D^S* (T - T.) 


+ R 4 (Tj a -Tf) -R 5 Tf, (i =i, 2 N) 

where 
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c. 

1 

T. 

1 
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C.. 
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total heat capacity of element i ( = p.v.c.) 

temperature of element i 

time 

conduction coupling term between elements i and j 



A.. 

11 

i.. 

ij 


Dl = a constant set equal to 0( lunar night) or 1 (lunar day) 

R = radiation coupling terms. They are, in the order listed in the 
equation, composed of the following quantities: 

(A Fa), (A Fa), (A F €€ <r ) , (A F'ee a ) 
tr i t s, a? m tmm 

and [ A ( 1 - F ) e c ] 
t m 

A^, = total and projected surface area, respectively, of element i 

F , F, F , F' , F = view factors with respect to reflected sunlight, direct in- 
r s, a m m * 

solation, shadow, and moon (with and without shadow) , respect- 
ively 

a,e = total solar absorptance and infrared emittance of vessel surface 

a = Stefan-Roltzmann's constant 

S = solar constant 

f = moon's total local albedo 

e = moon's infrared emittance 

m 
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The terms on the right-hand side of the equation represent, in the order listed, 
heat exchange with the following sources and sinks: adjacent elements, lunar 
reflected sunlight, insolation, shaded region of moon's surface, u n s haded re- 
gion of moon's surface, and space. Equation (i) is required for each element 
in order to compute the temperature as a function of time and is solved numeri- 
cally. 


Knowing each temperature, the heat transfer to the liquid per unit time 
is calculated by 

4 -£ c w<V V (2> 


where 


n = number of elements adjacent to liquid 
C. j= conduction coupling terms from element i to liquid [= 
= temperature of cryogenic liquid 


and the sum total of the energy transferred over any time increment is 
m 


k. .A. . 
l-l i-i 


i-i 


q = Yj q A t 
1 


(3) 


where 


period 


At = a single time increment 

m = number of time increments counting from beginning of storage 


The heat transfer to the liquid because of the insulation temperature is, 
unfortunately, not the total. Structural penetrations are always present and 
cause a degradation in the effectiveness of the insulation. This has been studied 
experimentally ( Table I) . In this analysis attempts to calculate the heat trans- 
fer from structural penetrations were abandoned, being too unrealistic with 
such an ill-defined design. Rather, this heat was accounted for by assuming that 
a given percent of the total heat absorbed by the liquid resulted from structural 
penetrations. This was made a parameter, and three values were used: 0, 20 
and 50 percent. Therefore, the total heat absorbed by the liquid is 

« = 4+ (-T^-) 4=4 (Trr) < 4 > 
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where 


x = fraction of heat absorbed resulting from penetrations 


TABLE I 

RESULTS OF EXPERIMENTAL INVESTIGATIONS OF 
HEAT FLOW DURING SPACE STORAGE THROUGH 
PENETRATIONS 


Reference 

Insulation 

Percent of Total 
Energy Transferred 
Through Structural 
Penetrations 

Remarks 





General Electric 
Company [i] 

Linde SI-62 
& NRC-2 

28-85 

3 spherical vessels 
(2, 4, & 8 ft. in di- 
ameter) 

1 cylindrical vessel 
(3 ft. in diameter by 
6 ft. in length) 

Boeing Company 
[2] 

Linde 

24 

Ellipsoidal 8 ft. 
diameter vessel 


and the sum total, since the beginning of the storage period, is 

Q = ^ QAt. (5) 


The evaporation rate m at any time is 


hi = Q/ L ( 6) 

v 


where 


L^ = latent heat of evaporation 
and the total amount evaporated, m, at any time, is 


m = £Q/ L V (?) 
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If no venting occurs, the temperature of the liquid will rise (Fig. 2) . During 
such a process the temperature at any time is given by 

T f =T i + ~" c" ' q (8) 


where 


T. , = initial and final temperature, respectively, during time increment A t 

C = heat capacity of liquid [ c ( m 0 -m) ] 

c = specific heat capacity of liquid 

m 0 = initial liquid mass 

m = liquid mass evaporated as given by eq. (7) 


If venting occurs, the temperature will remain constant sc long as the 
ullage pressure is constant ( Fig. 2) . The vent rate m* at any time is 


m* = 




(9) 


and the total amount vented m* up to time t is 



( 10 ) 


where 


P 



l 


m 

m 


density of cryogenic vapor and liquid, respectively 
evaporation rate 
total evaporated liquid 


VIEW FACTORS 
Diffuse Radiation 

The radiation exchange terms in eq. (i) require view factors for each 
surface element. The view factor for the direct impinging sunlight (Fig. 3) is 

F = cos a (11) 
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where 


a 


is angle between vectors N and S (Fig. 3) . 


Four view factors for diffuse lunar radiation are necessary: (1) reflec- 
ted sunlight, F , (2) shadow emission F (when such exists) , (3) emission 
r s, a 

from unshaded part of moon, F " , or (4) emission from total moon when no 

m 

shadow exists, F . These view factors are computed from the same expression 
m 

although modifications* are necessary, i. e. , F g ^ will be with respect to a 
much smaller area than F 


m 


_ cos y 
F = L - 


7 r R z 


The basic expression [3] is 

cot y cos' 


R 2 sin 


-» (-^r) 


( cot? y + A 2 

A 1 + R 2 / 


7T L 2 ( R+1 > 


R“ + 1 


(cot? y + 1 ) 2 


sin y 


it 


(R 2 -c ot?y) 2 . cos 1 
R 2 + 1 + 


cot? y + i\ 
1 + R 2 / 


i-i 

2 


( cot 2 y + 1 ) 




(12) 


where 


y = angle between normals to Vessel element and lunar surface 

R= ratio of radius of circular area representing the moon's surface to 

r 

height of vessel element above surface of moon ( = — — ) 

h 

This view factor applies for radiation exchange between a differential element 
and a finite circular area of radius r . The differential element is centered and 
above the circular area by a distance h. The elements on the vessel, because 


* Except for F which will be equal to F * for the diffuse case. 

r m 

* * For cases where R = °° this becomes, by L'Hospital rule, \ (1 + cosy ) or 

cos 2 (-f-) : if cot y > R, set cot y = R. For cases where y = 0, 


F = 


r 2 + h 2 
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they are very small compared to the circular area, are taken as the differential 
area. 


To get the view factor for the shaded part of the moon, the size of the 
shadow must be estimated. It will vary from a circle (with a diameter equal 
to the diameter of the storage vessel) , when the sun is directly overhead, to a 
rectangle of great length, when the sun is just rising or setting. This length 
( Fig. 4) is given by 


= h + 2r 
tan v 

The angle v is given by 


(13) 


v = sin 1 ( sin £ vi - sin 2 A ) 


(14) 


where 

l - time angle of sun ( equal to 0° when sun first rises above horizon, to 
90° when sun is directly overhead, etc . ) 

A = latitude location of vessel 

If the vessel is a cylinder with hemispherical ends, the height h of an element 
(above the moon's surface) located on the lower hemisphere, middle cylinder, 
and upper hemisphere are given by 

h = d + r + rN (15a) 

z 

h = d + r + Jh (15b) 

h = d + r + h + rN (15c) 

z 

respectively, where 

d = vertical distance from lowest point of vessel to lunar surface 
( = 1 for all calculations) 

r = radius of vessel 

h = height of middle 

if = vertical component of N ( N is the vector through center of element 
from center of vessel) 
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When the vessel is a sphere, eq. (15b) is not involved and when the vessel is 
flat ended, eqs. (15a) and (15c) are not involved. The ratio R , required for 

S 

calculating with eq. ( 12) the shadow view factor, F , is 

s 


R = 
s 


i 

2h 


( 16 ) 


Further modifications are required since the view factor calculated using eq. 
( 12) is for a circular area (which the shadow is not, except when the sun is 
directly overhead) . An approximation to the correct value can be made by 


F" 

s 



F 

s 


(17) 


where ^ 

A 0 = area of shadow if circular with radius — — 

z 

A 0 = area of shadow if rectangular with length l 

Elements on the vessel that are fully facing the shadow will be more strongly 
influenced than others. Those that are on the side away from the shadow, of 
course, will not be influenced at all. This is included in the view factor calcula- 
tions for a given element, by finding the cosine of the angle between the - S and 

N vectors and multiplying by F" or 

s 

F = cos a ' F' (18) 

s, a s 


where 


a * = angle between - S and N 

F * = view factor calculated using R 
s s 

The view factor for the unshaded part of the moon is found from 

F" = F - DiF (19) 

m m s, a 

where 


D 1 = constant having value of 1 during lunar day and 0 during lunar night 
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Nondiffuse Radiation 


The nondiffuse radiation is included in the calculations by modifying the 
view factors for infrared radiation and solar reflected radiation. Previously, 
these were computed from the same expression [eq. ( 12) ] . 

For diffuse radiation, the view factor is given by 

Fj_ 2 = / / - cos LOOS 8 dA.dAt (20) 

A, A, ’ ^ 

If the radiation is nondiffuse, an expression must be found which characterizes 
the radiation, as does Lambert's equation for diffuse radiation. For the infrared 
radiation, the expression of Pettit and Nicholson [4] is used: 

I = I 0 cos 2 / 3 e (21) 

and for the reflected sunlight, the expression of Hapke [5] is used (after being 
normalized to equal I 0 at its maximum value) : 

1 - * (cosT - :L0 * ) (b (a , g>) (», 


where 


I = radiant energy per unit time per unit area per unit solid angle 
I 0 = maximum value of I 


e = angle between normal to lunar surface and emitted or reflected 
beam of radiation 


i = angle between normal to lunar surface and incident beam of 
solar energy 

a = phase angle, angle between emitted or reflected beam and incident 
beam 


B ( a , g) 


tan a 

2g 


(1 _ e -g/tan a )(3 _ e -g/tana ) 



B (a ,g) 
g 



— r — ( r is the radius of a tube, A is the mean free path of light 

A 

ray through lunar surface, and g is assumed to be 0. 8. ) 
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Because the view factor is between a differential element A 2 and a finite area 
A l9 the integration over A 2 is eliminated. The view factor expression is now 
written in a general way as 



^ (e, <p ) cos Q dAi 
k-K r 2 


(23) 


and the function ^ ( e, cp ) is given by the expressions which follow. For diffuse 
radiation, cos e; for nondiffuse infrared radiation, cos 2 / 3 e ; and for nondiffuse 

reflected radiation, ( C0S 

\ cos 1 + cos e 

The k in the denominator is a normalizing factor equal to unity for the diffuse 
case and other than unity for the nondiffuse case. The view factors are calculated 
by numerical methods and substituted into the proper terms in eq. (1) . The 
results are presented later. 


)(- 


sing + ( 7T - a) cos 


IT 


— )(b (g,g)). 


COMPUTER PROGRAM 
General 

A computer program has been prepared to predict the thermal behavior 
of the vessel and its contents. The program is a combination of several smaller 
programs. The cross-sectional area and lengths for the conduction terms are 
calculated in a separate subroutine [6] . The view factors for nondiffuse radi- 
ation are computed in a separate program and read into the main program as in- 
put data. The view factors for diffuse radiation are calculated within the main 
program. The moon's temperature is computed from a Fourier series [7] 
which is a subroutine. * The shadow cast on the lunar surface was assumed to 
have a uniform temperature of 109 °K. Theoretical calculations have shown that 
no appreciable gradient exists in the region about the edge of the shadow and 
that a single constant temperature can be assigned to the shadow area [9] . The 
Heun method [10] with an automatic time step criteria [11] is used to solve the 
differential equations. Computer time for a run over one lunar cycle varies 
( 16 to 170 minutes) , but 80 percent of the runs took 20 minutes or less. A 
report is being prepared which gives the details of the program, i.e. , listing, 
flowchart, deck setup instructions, etc. [12]. 


* The Fourier series was fitted to Sinton's data [8], Variations in temperature 
with latitude are taken into account by multiplying by cos X where X is the lati- 
tude. 
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The computer program parameters are listed in Tables n and in. 

TABLE II 
INPUT DATA 


Parameter 

Identification of the number of elements in each layer of 
insulation ( Table VI) 

Insulation thickness 

Number of layers into which insulation is divided 
Radius of vessel 

Height of cylindrical middle, when present 

Designates shape; sphere, cylinder with flat ends, cylinder 
with hemispherical ends, respectively 

Latitude location of vessel 

Solar absorptanee of vessel surface 

Infrared emittance of vessel surface 

Insulation thermal conductivity normal to insulation 

Insulation thermal conductivity parallel to insulation 

Insulation specific heat 

Insulation density 

Cryogenic liquid specific heat 

Ciyogenic liquid density 

Vessel ullage gas density 

Albedo of moon 

Vent valve opening pressure 

Metric ( MKS) or English system of units 


Notation 

N 


L 

K 

r 

h 

0 , 1, -1 


A 


a 

e 

k * 

n 

V 

C* 



None 

None 


* These property data are temperature dependent. They are either read-in 
from a table (i. e. , Table IV) or expressions are used to give the values as 
a function of temperature. 
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TABLE IE 


OUTPUT DATA 


Parameter 

Notation 

Temperature of each element, of cryogenic liquid, and of 
moon 

T 

Cross-sectional area of elements, lengths between adjacent 
elements, volume of elements and other geometry data 
[6] 

A,i,V 

Insulation weight 

None 

View factors 

F 

Heat transferred to liquid through insulation and through 
penetrations ( supports) 

q 

Evaporation rate 

m 

Vent rate 

m* 

Percent vented 

None 

Initial liquid mass 

mo 

Liquid mass evaporated 

m 


Property Data 

The temperatures at the end of each time step were used as a basis for 
selecting the thermal properties for the succeeding computation. The thermal 
property data used for the multifoil insulation are presented in Tables IV and V. 
Table IV gives data that approximately coincide with the NRC-2 aluminized 
Mylar insulation. For NRC-2 insufficient data were available over the tempera- 
ture range of interest, therefore, assumptions* were necessary. The values 
that were taken from the literature (and they comprise only a few) are noted. 
Table V gives data taken from information supplied by the Linde Company on 
their super-insulation. No assumptions are involved here. 


* These assumptions involved extrapolations, interpolations, and assumptions 
about the slope of the curves. 
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The data used for the liquid hydrogen were taken from reference 13. All 
the data are for equilibrium liquid hydrogen. 

TABLE IV 

MULTIFOIL INSULATION DATA SIMILAR TO NRC-2 ALUMINIZED MYLAR 



(W/m -°K) 
0.082 
0.138 
0.148 [16] 
0.149 
0.150 
0.150 
0.150 


P [14] 
(kg/ m 3 ) 

41.648 

41.648 

41.648 

41.648 

41.648 

41 . 648 

41.648 


MULTIFOIL INSULATION DATA TAKEN FROM INFORMATION ON 
LINDE SUPER INSULATION * 


T 

( °K) 

c 

( J/kg - 0 K) 

kn 

(W/m - ° K) 

P 

(kg/m 3 ) 

16.7 

0.01255 x 10 3 

0.2767 x 10 -5 

112.1 

111.1 

0.3766 x 10 3 

0.1280 x 10“ 4 

112.1 

222.2 

0.7113 x 10 3 

0.4324 x 10~ 4 

112.1 

333.3 

0. 9623 x 10 s 

1.159 x 10" 4 

112.1 

444.4 

1. 172 x 10 3 

2.421 x 10~ 4 

112.1 

500.0 

1.297 x 10 3 

3.459 x 10" 4 

112.1 

833.3 

1. 883 x 10 3 

1.557 x 10~ 3 

112.1 


* Data obtained from the Linde Company. 
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Configurations 

The computer is programmed to perform calculations for three vessel 
shapes: a sphere, a cylinder with hemispherical ends, and a cylinder with flat 
ends. The insulation for each shape can be partitioned three different ways. 
Each way provides a different number of elements . The insulation can be fur- 
ther partitioned into layers. Three was the normal number of layers used for 
almost all calculations. This information is summarized in Table VI and is 
depicted in Figures 5 through 7. The information in Table VI is for one insu- 
lation layer; for three layers the number of elements will be tripled. 

TABLE VI 

IDENTIFYING NOTATION AND NUMBER OF ELEMENTS 
FOR ONE LAYER OF INSULATION 


N 

Sphere 

Cylinder with 
Hemispherical Ends 

N 

Cylinder with 
Flat Ends 


( Code = 0) 

( Code = -1) 


( Code = +1) . 

8 

8 

12 

8 

12 

32 

32 

40 

16 

24 

72 

72 

84 

24 

36 


Time Reference 


In this analysis, time (in hours) is given relative to the lunar midnight. 
The mean synodic period (time required for subsolar point to travel full 360° 
around the lunar surface) is 29.530589 earth days and is taken as one full lunar 
cycle. Calculations requiring time dependent angular measurements are rel- 
ative to the time angle £ defined by the following expression 

£ J -3SS- )t 

? \( 29. 530589) (24) / 

£ = 0. 507948t (24) 

where t is in hours. 
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RESULTS 


General 

Table VII shows the numerical values used to make the calculations. 

The values apply to most, though not all, runs (the exceptions are noted) . Run 
5 is typical of the runs for which the values apply. Table VIII summarizes the 
computer results. Table VIII (reading from left to right) presents the following: 
the run number, identification for the number of elements in each layer of insu- 
lation (Table VI) , the code that designates the tank shape (Table VI) , the radius 
of the vessel, the lunar latitude location of the vessel, the percent of the liquid 
vented for an idealized case (no heat enters through penetrations) the ratio of 
percent vented to percent vented for Run 5 (used as a standard of comparison), 
and pertinent remarks that distinguish one run from the other. The initial 
mass of liquid hydrogen was obtained by using tank dimensions which correspond 
roughly with those of payloads compatible with the Lunar Module. This approach 
was taken by A. D. Little [9] and seems reasonable. However, one other size 
has been investigated ( Table VIII, Run 2) . 

Temperatures 

Figure 8 shows the temperatures for Run 5 for the elements indicated. 

The location of the elements are shown in Figure 5a, pertinent information 
about the run is given in Table VIII. * ** Elements 5 and 6 are in the third layer. 
The temperature characteristics are reasonable ( and typical of the other runs) 
when the location of each element is considered. Figure 9 shows the temperature 
for a case where calculations were based on a very thin layer of insulation. 

This was done to allow the temperatures to achieve steady state. These temper- 
atures were used, with the regular insulation thickness (10. 16 cm; 4 inches) , to 
calculate the venting. Although the temperature difference between Run 5 (out- 
side layer) and Run 10 is negligible, some difference is present in the heat 
transfer (Fig. 10) . The shift in the curves is caused by the differences in in- 
sulation heat capacity ( much higher for Run 5) . If the curves peaked at the 


* An estimate of the venting for cases where structural penetrations are present 
may be made by straight multiplication of the values in the table. For example, 
if 50 percent Q f the total flux is assumed to be due to penetrations, the table 
value for percent vented would be doubled. 

** Information about all figures can be obtained by referring to Table VIII. 
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TABLE VII 


L 

m 0 


r 


e 

m 

f 

S 

€ 

a 

T 

P 


t 


A 


t 


<j 


NUMERICAL VALUES USED IN COMPUTATIONS 


Insulation weight 1 

223.3 kg (492.3 lbs) 

Insulation thickness 2 

10.16 cm(4. 0 inches) 

Initial propellant 3 
weight 

2187.9 kg (4823.4 lbs) 

Vessel radius 4 

2.1m (6.88 ft) 

Vessel volume 5 

38.79 m 3 (1369. 86 ft 3 ) 

Moon’s infrared emittance 

1.00 

Albedo of moon 

0.07 

Solar constant 

1395 W/m 2 

Vessel infrared 6 
emittance 

0. 90 

Vessel solar absorptance 6 

0.20 

Initial liquid hydrogen 
temperature 

20. 4°K 

Initial ullage pressure 

1 atm 

Initial ullage 

10 % by volume 

Radius of lunar surface 

520 m (1706 ft) 

Projected surface area 7 
of an element on vessel 

4.41 m 2 (47.47 ft 2 ) 

Surface area of an element 7 
on vessel 

6. 93 m 2 (74. 59 ft 2 ) 

Stefan-Boltzmann’s constant 

5. 673 x 10 -8 W/m 2 - °K 4 


*A11 weights are with respect to earth gravity. Run 2, 53.13 kg (117.1 
lbs); Run 21, 601 . 1 kg ( 1325. 2 lbs) ; Run 10, 4.61 kg (10.16 lbs) ; Run 
3, 269. 8 kg ( 594. 8 lbs) ; Run 4, 233. 7 kg ( 515. 2 lbs) . 


2 Run 22, 3. 534 cm (1.39 in) ; Run 10, 0. 2 cm ( 0. 07874 in) 

3 Run 2, 273. 5 kg (603. 0 lbs) 

4 Run 2, 1.05 m (3.44 ft) ; Run 4, 1. 68 m ( 5. 51 ft) 

5 Run 2, 4. 85 m 3 (171.24 ft 3 ) 

6 Run 14, a = e = 1. 0 

7 Does not apply to Runs 2, 3, 4, and 7. Run 7, A = 0. 74 m 2 (0. 80 ft 2 ) , 
A = 0. 77 m 2 (0. 83 ft 2 ) t 


17 



TABLE VIE 

SUMMARY OF CALCULATIONS 


Run 

No. 

N 

Shape 

Code 

r 

X 

Percent 
Vented 
(1 lunation) 

Percent 

Vented 

Ratio 

Remarks 

( Deviations from Run 5) 

1 

8 

0 

2.1 

30 

5.53 

0.96 


2 

8 

0 

1.05 

30 

11.39 

2. 06* 


3 

8 

+1 

2.1 

30 

7.35 

1. 33* 

h = 2. 8 

4 

8 

-1 

1.68 

30 

6.17 

1. 12* 

h =2.13507 

5 

8 

0 

2.1 

0 

5.79 

1.00 


7 

72 

0 

2.1 

30 

6.20 

1.12* 

216 elements 

8 

8 

0 

2.1 

0 

22. 09 

3. 82 

Property data constant 
at 300 °K 

9 

8 

0 

2.1 

0 

5. 81 

1.00 

ii 

o 

10 

8 

0 

2.1 

0 

4. 55 

0.79 

Thin insulation, k = 0 

11 

8 

0 

2.1 

0 

6.64 

1.15 

Shadow view factor is 
zero 

14 

8 

0 

2.1 

0 

8.51 

1.47 

a. = e = 1 

15 

8 

0 

2.1 

0 

9.46 

1.63 

Property data constant 
at 200 °K 

19 

8 

0 

2.1 

0 

5. 89 

1.02 

cos 2 / 3 e function 

20 

8 

0 

2.1 

0 

6.58 

0 . 99** 

Hapke function, cos 2 / 3 e 
function, shadow view 
factor is zero 

21 

8 

0 

2.1 

0 

0.37 

0.06 ; 

Linde insulation with 
thickness 10. 16 cm 

22 

8 

0 

2.1 

0 

1.45 

0.25 

Linde insulation with 
thickness 3. 534 cm 

23 

8 

0 

2.1 

0 

6.64 

1 . 00** 

Hapke function, shadow 
view factor is zero 

24 

8 

0 

2.1 

0 

6.53 

0.98** 

cos 1 / 9 e function, shadow 
view factor is zero 

25 

8 

0 

2.1 

0 

6.58 

0 . 99** 

cos 2 / 3 e function, shadow 
view factor is zero 


* Compared with Run i 
** Compared with Run 11 
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same thermal energy value, the sums, if taken over a lunation, would be about 
the same. For this reason some investigators [9, 17] have used the steady 
state method to obtain satisfactory answers for this type of problem. In this 
investigation the thermal energy absorbed by the liquid hydrogen was generally 
higher for Run 5 than for Run 10 (Fig. 10) . This is attributed to three factors: 
the insulation thickness, which was less for Run 5 since the overall insulation 
thickness is divided into 3 layers ( the flux is computed using the thickness of 
the innermost layer ) ; the temperature difference, which is less for Run 5 
since the temperature of the innermost layer is used; the thermal conductivity 
of the insulation, which is lower for Run 5 since the temperature is lower and 
k^, being temperature-dependent, decreases. The first factor increases the 

flux and the last two cause a decrease. The net result during most of the lunar 
cycle (the lunar morning excepted) is a greater thermal energy transfer for 
Run 5. 


Thermal Conductivity 

The percent vented ratio ( Table VIII) clearly demonstrates the influence 
of a number of factors on the calculated results. The factor exerting the greatest 
influence is the thermal conductivity and is aptly demonstrated by a comparison 
Of Runs 21 and 5. The insulation data shown in Table V were used in Run 21; 
the data in Table IV were used in Run 5. The values of k^ for Run 21 are rough- 
ly an order of magnitude smaller than for Run 5, The insulation weight increases 
by about a factor of 3 ( Table VII) . Run 22 was made keeping the insulation weight 
equal to that of Run 5 (varying the thickness) and the results still indicate the 
overriding influence of k . Run 9 shows the negligible influence on the calculated 
results when the thermafconductivity in the lateral direction (k^) is zero. 

Thermal properties, if assumed to be temperature independent, can lead 
to significant errors. This is shown by Runs 8 and 15 in Table VIII. 

Shadow Cast on Moon' s Surface 

The effect of the shadow is not negligible. Disregarding the shadow 
(Run 11) gives a 15 percent increase in the percent vented over the value in 
Run 5 where it is included. The appreciable temperature difference between 
Run 5 and Run 11 is further evidence of its significance (Fig. 11) . 

View Factors 

The view factor calculations for diffuse radiation are presented in Figures 
12 through 15. Equation (12) was used to compute the view factors for infrared 
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and reflected radiation from the lunar surface. In Figure 12 the shadow view 
factor is hardly noticeable, having a maximum at 0.73 lunation (lunar noon 
occurs at 0.50 lunation) . The location of this element (Fig. 5a) readily explains 
why. Figure 13 gives the same information for another element. 

Figures 14 and 15 give the view factors for elements located on a cylindri- 
cal vessel with flat ends located at a latitude 30 degrees above the lunar equator 
(previously the vessel was a sphere located on the equator) . In Figure 14 the 
element is facing the surface and in Figure 15 the element is vertical to the 
surface. 

The view factor calculations for nondiffuse radiation are presented in 
Figures 16 and 17. For comparison, the same elements used in Figures 12 and 
13 were plotted. Equation (23) was used to compute the view factors for infrared 
and reflected radiation from the lunar surface using the cos 2 / 3 e and Hapke ex- 
pressions, respectively. The shadow view factor was excluded from these cal- 
culations so the effect of reflected radiation could be determined when at its 
maximum (the shadow will attenuate radiation that has a large backscatter com- 
ponent) . Figures 16 and 17 show that the view factors for reflected sunlight 
are very strongly dependent on the position of the sun (sunset occurs at 0. 75 
lunation) . The view factors for infrared radiation are only slightly different 
from the diffuse case, and because of this, cos 9 e was arbitrarily picked 
for comparative purposes to make additional calculations. The results for all 
three cases (cos e, cos 2 / 3 e, and cos V 9 e) are given in Table IX. The elements 
are located on a vessel as in Run 5. 

TABLE IX 

VIEW FACTORS FOR NONDIFFUSE 
INFRARED RADIATION FROM LUNAR SURFACE 


Element 

cos e 

cos 2 / 3 e 

cos 1 / 9 e 

1,1,1 

0.18 

0.21 

0.27 

2,4,1 

0. 82 

0.79 

0.74 


Notice in Table IX that for element 1,1,1 (upper part of vessel) the 
values increase from case to case, while for element 2, 4, 1 (lower part of 
vessel) the values decrease. Also the sum of the two elements in each case is 
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approximately unity. This is logical because the total radiation striking the 
vessel does not change (as long as the temperature of the lunar surface at the 
location of the vessel is constant) , only the distribution of the radiation. 

Nondiffu se Radiation 

The influence of nondiffuse radiation from the lunar surface is included 
in the calculations, through the nondiffuse view factors. This effect on the tem- 
perature of the vessel elements, though small, is present (Figs. 18 and 19) . 

Run 11 is diffuse infrared radiation given by Lambert's cos e expression and 
Run 25 is nondiffuse infrared radiation given by Pettit and Nicholson's cos 2 / 3 e 
expression. 

Notice that the temperature for element 4 is lower in Run 25 than in 11 
(Fig. 19) . The reverse is true for element 1 (Fig. 18) . This follows from the 
view factors in Table IX. Because cos 2 / 3 e is such a slight departure from the 
diffuse case, as evidenced by the results, Run 24 was made using cos 1 / 9 e. 

This expression was arbitrarily picked for comparative purposes and has no 
theoretical or experimental basis. 

The reflected radiation in Runs 11 and 25 was assumed to be diffuse. 
Hapke's expression [eq. (22)] was used for nondiffuse reflected sunlight. Al- 
though the view factors are strongly affected, the temperatures are not because 
of the low albedo (0. 07) of the moon. This is shown in Figures 20 and 21. In 
Run 20 the cos 2 / 3 e and Hapke’s expression have been combined. Runs 20 and 
11 are compared in Figure 22. 

Shadow effects were not included in Runs 11, 20, 23, 24, and 25. Runs 
5 and 19 contain shadow effects. Figure 23 shows that for an element facing 
the lunar surface the nondiffuse effect is washed-out by the shadow ( except 
during the lunar morning) . This is not so for elements facing away from the 
lunar surface as shown in Figure 24. 

Table VIII shows that there is no influence of nondiffuse radiation on the 
venting. This is because the vessel receives the same total amount of energy, 
since it receives radiation from all directions, regardless of the spatial distri- 
bution of the radiation from the lunar surface. 

Shape and Size of Vessel 

Venting is influenced by the shape of the vessel as shown by Runs 1, 3, 
and 4 (Table VIII) . All vessels have the same volume. Venting is also a strong 
function of size, increasing as the vessel becomes smaller (Run 2) . 
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Number of Isothermal Elements 


The accuracy of the calculations can be increased at the expense of over- 
all complexity by dividing the insulation into a greater number of elements. This 
has been done in Run 7. The insulation is divided into 216 elements (most cal- 
culations are for 24 elements ) and the 12 percent difference between Runs 1 
and 7 represents the accuracy sacrificed when the more simple case (Run 1) is 
used. Of interest in these two cases is a temperature comparison (Fig. 25) for 
elements of similar location (Figs. 5a and 5c) but different size. 

Absorptance and Emittance 

The a/ e values of the outer surface have a profound effect on the venting 
and in this analysis are assumed to provide, along with the insulation, the only 
means of thermal control. The accumulation of materials from the moon’s sur- 
face onto the vessel surface could alter the a / e ratio. For this reason, Run 
14 was made where a = € = 1 was used. The results are given in Table VIII. 

CONCLUSIONS 


The following factors were found to have a significant influence on the 
thermal characteristics of a storage vessel in a lunar simulated environment: 
thermal conductivity normal to insulation, constant thermal properties, size 
and shape of the storage vessel, a/ e ratio, and the shadow cast by the vessel. 

Nondiffuse infrared radiation (on temperature only) , the number of iso- 
thermal elements on the vessel surface, and steady-state temperatures (on 
venting only) , were found to have slight, though not necessarily negligible in- 
fluence on the vessel. 

There was no appreciable influence resulting from the lateral thermal 
conductivity of insulation and reflected sunlight (neither diffuse nor nondiffuse) . 


22 




FIGURE 1. THREE BASIC SHAPES 
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5 6 7 8 9 10 II 12 

ENTROPY cal/gm°K 


1. No venting (1-2) 

2. Venting begins and pressure and temperature remains constant ( 2-3) 

3. Venting occurs during entire storage period (1-4) 


FIGURE 2. TEMPERATURE - ENTROPY DIAGRAM FOR VENTING 


23 







z 


z 




(a) 


(c) N=72-NUMBER OF ELEMENTS =84 


FIGURE 6. CYLINDER WITH HEMISPHERICAL ENDS (CODE NUMBER = -1) 
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(o) N-8- NUMBER OF ELEMENTS = 12 (b) N= 16- NUMBER OF ELEMENTS 

Z 



(c ) N * 24— NUMBER OF ELEMENTS = 36 


FIGURE 7. CYLINDER WITH FLAT ENDS ( CODE NUMBER = +1) 






TEMPERATURE 



TIME FRACTION OF ONE LUNATION 


1 . Zero time corresponds to the lunar midnight. 

2. One lunation is equal to 708. 7 hours or 29. 5 earth days. 

3. Element location is given in Figure 5a. 
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FIGURE 8. SPHERICAL TANK ELEMENT TEMPERATURE 
VARIATIONS (RUN 5) 







TEMPERATURE °K 


700 



TIME FRACTION OF ONE LUNATION 


1. Zero time corresponds to the lunar midnight. 

2. One lunation is equal to 708. 7 hours or 29. 5 earth days. 

3. Element location is given in Figure 5a. 


FIGURE 9. ELEMENT TEMPERATURE VARIATIONS BASED ON 
A THIN LAYER OF INSULATION (RUN 10) 
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TIME FRACTION OF ONE LUNATION 


FIGURE 12. VIEW FACTORS FOR ELEMENT 1, 1, 1 (RUN 5) 
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FIGURE 14. VIEW FACTORS 



OF ONE LUNATION 


, 4, 1 ( RUN 3) 







VIEW FACTOR 



TIME FRACTION OF ONE LUNATION 


FIGURE 17. VIEW FACTORS FOR NONDIFFUSE RADIATION (RUN 20) 

ELEMENT 2, 4, 1 
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TIME FRACTION OF ONE LUNATION 

FIGURE 1 8. TEMPERATURE COMPARISON BETWEEN DIFFUSE COS e 
( RUN ii), NONDIFFUSE COS 2 / 3 e (RUN 25) AND COS 1 / 9 e (RUN 24) 
INFRARED RADIATION. ELEMENT 1,1,1 





TEMPERATURE 



TIME FRACTION OF ONE LUNATION 


FIGURE 20. TEMPERATURE COMPARISON BETWEEN DIFFUSE 
(RUN 11) AND NONDIFFUSE (RUN 23) REFLECTED RADIATION. 

ELEMENT 1,1,1 


40 





TIME FRACTION OF ONE LUNATION 


FIGURE 22. TEMPERATURE COMPARISON BETWEEN DIFFUSE (RUN 11) 
AND NONDIFFUSE (RUN 20) FOR COMBINED INFRARED AND 
REFLECTED RADIATION. ELEMENT 1,1,1 




TEMPI 





FIGURE 24. TEMPERATURE COMPARISON BETWEEN DIFFUSE (RUN 5) 
AND NONDIFFUSE (RUN 19) INFRARED RADIATION. 
ELEMENT i, 1, 1 
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TIME FRACTION OF ONE LUNATION 


FIGURE 25. TEMPERATURE COMPARISON BETWEEN ELEMENTS OF 
DIFFERENT SIZE BUT SIMILAR LOCATION ON THE 
STORAGE VESSEL 
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